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ABSTRACT 

NGC 663 is a young, moderately rich open cluster, known to contain one of the largest frac- 
tions of Be stars among all galactic clusters. In this work we present CCD uvby/3 photom- 
etry for stars in its central area. We have used these data to obtain the main cluster phys- 
ical parameters. We find that the reddening is highly variable, with values ranging from 
E(b -y)^ 0.639 + 0.032 in the centi-al part to E{b - y) = 0.555 + 0.038 in the south-east. The 
distance modulus is found to be 1 1.6 + 0.1 mag. (2.1 Kpc), and the age log t - 7.4 + 0.1 years 
(25;^^ Myr). The age obtained is consistent with the interpretation of the Be phenomenon as 
an evolutionary effect. 

Key words: techniques: photometric - stars: emission-line. Be - Hertzsprung-Russell (HR) 
diagram - open clusters and associations: individual: NGC 663. 



1 INTRODUCTION 

The precise determination of the main physical parameters of 
galactic open clusters plays a central role in the study of the stellar 
structure and evolution. With accurate photometric data, and once 
the external variables such as reddening are corrected for, the clus- 
ter distances, ages and chemical abundances can be inferred from 
the study of the photometric colour-magnitude and colour-colour 
diagrams. 

This work is part of a programme aimed to produce accurate 
and homogeneous dating for a sample of young galactic clusters by 
means of CCD Stromgren uvby/3 photometry. The use of isochrone 
fitting in the Vq - cq plane of the iivby photometric system is an ade- 
quate tool to obtain accurate cluster ages. The range of variation of 
the Co index along the B-type sequence amounts to more than one 
mag., being significantly larger than most commonly used photo- 
metric colours. Moreover, the cg index is less affected by reddening, 
and allows an efficient segregation of emission-line stars. 

NGC 663 is a moderately rich open cluster in Cassiopeia, lo- 
cated in the Perseus galactic arm. It is thought to be the core of the 
Cas OB8 association. A remarkable characteristic is that it con- 
tains one of the largest fractions of Be stars among all galactic 
clusters. More than 30 Be stars have been detected so far among 
its B star sequence, and at least 30 per cent of all early B clus- 
ter members are known to be Be stars i lSanduleakllT979lll990l: 
[Pieulski. KoDacki & Kolaczkowski 2001gj). 

The accurate age determination of clusters with a high Be star 
content is a key element of the on-going discussion about the evo- 
lutionary status of Be stars. Fabreaat & Torreion (2000) proposed 
that the Be phenomenon is an evolutionary effect which appears at 
the end of the B stars main sequence lifetime. The main argument 
in favour of this hypothesis is their preliminary finding, based on an 
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inhomogeneous set of age determinations from the literature, that 
clusters with high Be star content have ages in the interval 13-25 
Myr, while no - or very few - Be stars are found in clusters younger 
than 10 Myr. Hence, the accurate determination of the age of NGC 
663 is an important issue in this discussion. 

NGC 663 has been the subject of many photometric 
studies in the past. Photoel ec tric UBV photome tr y wa s ob- 
tained by iHoaeetalJ <196ll). iMcCuskev & Hou3 il964 and 
Ivan den Bergl^^d^Roual^T^) . Stromg ren uvby and Crawford- 
Barnes H/? photometry was presented by iTapia et SI il99 l') and 
l^abregat et al. ( 1996). The former also obtained infrared JHK pho- 
tometry. Recent studies based on CCD photomet ry in the lohn- 
son U BV and Cousins RI systems include those of'Phel ps & lane j 
(1991) and Pigulski. Kooacki & Kolaczkowski (2001^. Photo- 
metric Ha surveys aimed to detect new Be stars have be en con- 
ducted by .Capilla. Fabregat & Baines, 12000.) and .Pigulski et alJ 
j200lj) . 

The cluster physical parameters determined in the above 
studies show an important dispersion. A reason for this could 
be the strongly variable reddening across the cluster area (e.g. 
|yan den Bergh & de Roux 1978) which makes the determination 
of the intrinsic magnitudes and colours more difficult. The 
publi shed distance mod uli range from 11.4 mag. (1.9 Kpc, 
McC uskev & Houkll964 to 12.25 mag. (2.8 K pc.lPhelps & lane j 
11994) . The determined ages range from 9 Myr. fl^pi^^l^^^]) 
to 20-25 Myr. ( Pigulski et al. 2001a) . 

In this work we present CCD uvby/3 photometry for stars in the 
central area of NGC 663 and a precise determination of the cluster 
main physical parameters. A detailed study of the cluster Be star 
population will be presented in a forthcoming paper. 
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Table 1. Exposure times with the different filters used, in seconds. 



Filter 


short 


medium 


large 


y 


10 


50 


200 


b 


12 


60 


240 


V 


35 


175 


700 


u 


120 


600 


2x1200 




12 


60 


240 




30 


150 


600 



Table 2. List of observations 



JD date airmass exposure 



2451 138 


20-11-98 


1.26 


short 


2451 138 


20-11-98 


1.40 


long 


2451 139 


21-11-98 


1.28 


short 


2451 139 


21-11-98 


1.46 


long 


2451 140 


22-11-98 


1.15 


short 


2451 140 


22-11-98 


1.17 


long 


2451 140 


22-11-98 


1.45 


medium 


2451 140 


22-11-98 


1.56 


medium 



2 THE DATA 

CCD photometry of NGC 663 was obtained on the nights 20 to 
22 November 1998 at the 1.52-m. telescope of the Observato- 
rio Astronomico Nacional, located at the Calar Alto Observatory 
(Almeria, Spain). The chip employed was the Tektronics TK 1024 
AB, with a size of 1024 x 1024 pixels. The 0'.'4 unbinned pixels pro- 
vide a field size of 6'.9x6'.9. Two slightly overlapping fields near 
the cluster centre were observed to cover the cluster central region. 

Observations were done through the four Stromgren uvby and 
Crawford narrow and wide Hy8 filters, every field being sequentially 
measured through the six filters. Three different exposure times 
were used with each filter, in order to ensure a wide range of stellar 
magnitudes. Exposure times in each filter were selected so that a 
B type star produces approximately equal count rates through all 
filters. Employed exposure times are presented in TableQ The list 
of all observations is presented in Table|2| 

During the same observing run we also obtained photometry 
of th e clusters h & x Persei, which has been published elsewhere 
iCap illa & Fabregat 2002, hereafter referred to as CF02). A de- 
tailed description of the image processing procedures, atmospheric 
extinction calculation and standard uvbyfi transformations is given 
in CF02. The accuracy of the standard photometry is 0.023, 0.014, 
0.017, 0.017 and 0.021 mag. in V, {b-y), m\, ci and j8 respectively. 
In CF02 it is also shown that our photometry is well tied to the 
standard system and free of systematic effects. 

Astronomical coordinates for all observed stars were derived 
from the instrumental ^ixel coord inates by using 26 stars from 
the Tycho 2 Catalogue rH0s et al."2000') included in the observed 
fields. Transformation equations were computed by means of the 
Starlink program ASTROM I Wallace 1998). The final astrometric 
accuracy, measured as the RMS of the mean catalogue minus trans- 
formed values for the stars used in the transformation, is better than 
0'.' 1 both in Right Ascension and Declination. 

Equatorial coordinates and mean photometric magnitudes, 
colours and indices for stars in the central region of NGC 663 are 
presented in Table |3] We have adopted the cluster star numbering 
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Figure 1. Colour-magnitude diagram for NGC 663. Filled and open circles 
denote stars considered as cluster members and nonmembers respectively 




^1 

Figure 2.V - c\ photometric diagrams. Symbols as in Fig.lTI 

system from the WEBDA database 1 fMermilliod"200(f). For clus- 
ter numbers lower than 635, WEBDA numbers are coincident with 
IWallenguislI i 19291) numbers, which we will refer to as 'W' here- 
after. 



3 REDDENING, INTRINSIC COLOURS AND DISTANCE 

The colour-magnitude diagram of all observed stars is presented in 
Fig.Q and the photometric V - ci diagram in Fig.|2| To obtain the 
intrinsic colours we first classified the stars as belonging to the early 
(earlier than AO), intermediate (AO- A3) or late (A3 onward) groups 
defined by Stromgren 1 1966). The cl assification was perfo rmed by 
means of the algorithm described bv lFigueras et alJil99lh . 

Cluster membership was assigned from the position of each 
star in the different uvbyfi photometric diagrams. Stars considered 
as members are marked in the last column of Table|3] and are rep- 
resented with a different symbol (filled circle) in Figs.Qand|5| The 
reliability of the photometric membership criteria was checked in 
CF02, where we conclude that our cluster photometric sequences 
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Table 3. Coordinates and photometric data for stars in tlie NGC 663 area. N is the number of uvhy and Hfi individual measurements. Stars considered as 
cluster members are denoted with 'm' in the last column. Only a few entries are displayed. The complete table is available in electronic form. 



Star 


a()2mO) 


(5(J2000) 


V 


b-y 


m\ 


c\ 


(TV 


O-b-y 


0"mi 


o-,, 


N 






N 


M 


1 


01 46 33.806 


+61 13 34.04 


14.653 


0.545 


-0.144 


0.819 


0.016 


0.017 


0.045 


0.065 


5 


2.767 


0.060 


3 


m 


2 


01 46 35.604 


+61 13 39.08 


12.258 


0.572 


-0.196 


0.230 


0.010 


0.007 


0.020 


0.021 


6 


2.473 


0.017 


3 


m 


3 


01 46 39.876 


+61 14 02.51 


13.029 


0.508 


-0.124 


0.538 


0.003 


0.009 


0.016 


0.022 


6 


2.595 


0.042 


3 


m 


4 


01 46 38.997 


+61 14 06.06 


11.025 


0.490 


-0.166 


0.131 


0.002 


0.016 


0.036 


0.034 


4 


2.558 


0.026 


3 


m 


5 


01 46 37.352 


+61 14 17.48 


14.288 


0.561 


-0.150 


0.663 


0.015 


0.016 


0.032 


0.036 


6 


2.708 


0.009 


3 


m 



Table 4. Mean reddening values for the three regions of the NGC 663 nu- 
cleus defined in Fig.l3l 



Region 


E{b-y) 


stars 


A 


0.590 ± 0.026 


100 


B 


0.639 ± 0.032 


40 


C 


0.555 ± 0.038 


49 



are not significantly contaminated with the inclusion of nonmem- 
ber field stars, and, conversely, very few, if any, actual members 
have been excluded from them. 

All member stars belong to the early group, indicating that 
our magnitude limited photometry does not reach beyond the late 
B spectral types. Reddening values and intrinsic colours and in- 
dices were henc e obtained by means of the procedure described by 
ICrawforj il978") . W e have used the standard (b - y)o - co relation 
given m Table VI of iPerrv. Olsen & Crawford h987h . and the fol- 
lowing reddening relations: 

Av = 4.3E(b - y) 
E{ct) = 0.2E(b - y) 
£(mi) = -Q.32E(b-y) 

Known supergiant and Be stars have not been included in the 
computation of the interstellar reddening. B type supergiants do not 
follow the standard (h - >')o - cq relation, while Be stars present an 
additional reddening contribution of circumstellar origin. 

We obtained a mean reddening value of E(b - y) = 0.591 ± 
0.042, from 189 stars. The large value of the standard deviation in- 
dicates the presence of variable reddening across the cluster area. 
In Fig.|3|we have represented the reddening values for individual 
stars as a function of their position. The cluster centre is the most 
heavily reddened region, while in the south-east part the extinc- 
tio n is sign ificantly lower than the cluster average, as already noted 
bv lFabregat et al. ( 1996). In Fig.|3|we have divided the cluster nu- 
cleus area into three regions of different reddening. Mean redden- 
ing values in these regions are presented in Table |4| Within each 
region, and in particular in regions B and C, the standard devia- 
tion of the mean reddening is still large, indicating that we cannot 
consider constant extinction within the region. Instead, the extinc- 
tion appears to be clumpy and irregularly variable with high spatial 
frequency. 

In Fig. we present the intrinsic Vo - cg diagram for B stars. 
Each star has been dereddened on the basis of its position within 
Fig. and using the reddening values in Tabled To obtain the 
distance we have fitt ed to the Vn - Cn d iagram the ZAMS as pre- 
sented in Table VI of IPerrv et sii Jl987h . We found the best fitting 
at a distance modulus of 1 1.6 mag. To estimate the error of this de- 
termination we have also represented in Fig.|3] as dotted lines, the 
ZAMS shifted by distance moduli of 1 1 .4 and 11.8 respectively. We 
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Figure 3. Reddening spatial distribution for B stars in NGC 663. Positions 
are relative to star W22 (0l'"46'"3a:24, +61° 12'59'.'6, J2000) 



find that the 11.6 value produces a distinctly better fitting than the 
two latter ones, and hence we give the value of 1 1.6 ± 0.1 mag. as 
the distance modulus of NGC 663. This value is in good agreement 
with the recent determinations based on CCD photometry. 



4 CLUSTER AGE 

Age determination has been done by isochrone fitting to the upper 
main sequence. In the uvhy system, for early-type stars the (b - y)o 
colour and the cq index are temperature indicators, and hence both 
Vq - (b - y)o and Vo - Cg planes are observational HR diagrams. 
Following the discussion in Fabresat & Torreion (2000) we con- 
sider the isochrone fitting to the Vq - Cq diagram as more precise 
and reliable, for the following reasons: i./ the range of variation 
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Figure 4. ZAMS fitting to the NGC 663 B star sequence. The soHd fine 
represents the ZAMS at a distance modulus of 11.6 mag. Dotted lines are 
for distance moduli of 1 1 .4 and 11.8 respectively. 
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Figure 5. Compaiison between the B star sequences of NGC 663 and h & 
X Persei . 



of the Co index along the B-type sequence is more than ten times 
larger than the range of variation of (b - y)o, providing much better 
discrimination between isochrones of similar ages; ii./ the co index 
is less affected, by a factor of 5, by interstellar reddening; iii./ the 
^0 - Co plane allows an efficient segregation of emission line stars. 

As a first approach, in Fig.|5|we compare the B star sequence 
of NGC 663 in the Vo - cq plane with that of h and x Persei. Pho- 
tometric data for h and x Per have been obtained from CF02. The 
turn-off of both sequences are clearly seen, and from them it is 
apparent that NGC 663 is significantly older than h and x Persei. 
Five NGC 663 stars are distinctly outside the cluster sequence. Star 
W4 lies at the left of the cluster turn-off, close to the ZAMS, and 
it is hence a cluster blue straggler. Its spectral type, derived from 
its C| index, is BIV. This star is also a 13 Cephei pulsating variable 
jPiet rzvfiski' 1 9971 IPigulski et a l. 2001bV Its position in the pho- 
tometric diagram is close to W161, another ma in sequence blue 
straggler already noted by iFabreeat et alj il996l) . but outside the 
field observed in this study. 
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Figure 6. Isochrone fitting to the NGC 663 B stai' sequence, in the Vq - co 
plane. 



The other four stars are in an horizontal sequence with My 
between -6 and -5 magnitudes. They are stars W40 (BD +60 
343), W44 (BD -^60 339), W54 (BD +60 333) and W86 (BD +60 
331). All of them are supergiants. They will not be considered for 
isochrone fitting in the My - co plane since the cq index defines 
different effective temperature scales for main sequence and super- 
giant stars. In CF02 there is a detailed discussion on this subject. 

In Fig. |6| we present the Vq - Cq sequence together with 
isochrones with ages of logf = 7.3, 7.4 and 7.5 years. The 
isochrones have been computed with the evolutionary models of 
ISchallereTal. 1 199 2), and transforme d to the observa tional plane 
by means of the relations obtained bv lTorreionI ^993)- Practically 
all stars at the cluster turn-off lie between the log? = 1 A and 7.5 
isochrones. Most stars are located on the 7.5 isochrone, but a sig- 
nificant number, 5 stars, are at the left, making the 7.4 isochrone 
better for defining the lower envelope of the stellar sequence. 

In order to use the supergiant stars to further restrain the age 
determination, in Fig. we present the Vq - {b - j)o photometric 
diagram, {b-y) is a measure of the Paschen continuum slope, which 
is correlated with effective temperature for stars of all luminosity 
classes. In this plane the four supergiants are closer to the log t = 
7.4 years isochrone than to the 7.5 one. Their positions even suggest 
a still younger age. However, for the reasons discussed above, the 
error in the position of a star in the Vq - (b- y) plane is likely to be 
affected by much bigger errors, and hence we give a higher weight 
to the age determination through the Vq - Co plane. 

From all this discussion we finally adopt an age of log t = 
7.4 ± 0. 1 years (25^j Myr). The quoted error is higher than the one 
we estimated in our age determination of the h and ^ Persei clusters 
in CF02. The reason for this is the more variable and clumpy nature 
of the interstellar reddening in the NGC 663 field, which makes the 
reddening determination more imprecise and the sequences in the 
intrinsic photometric planes more disperse. 



5 DISCUSSION 

iFabregat & Torre ion* proposed that the Be phenomenon is 

an evolutionary effect, appearing at the end of the main sequence 
lifetime of a Be star. They based their conclusion in the analysis of 
the frequency of Be stars in clusters of different ages. They found 
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Figure 7. Va-{b - y)o plane for B type stars in NGC 663. 



that in clusters younger than 10 Myr Be stars are lacking or very 
scarce. The maximum Be star frequency occurs in clusters in the 
13-25 Myr age interval. The main sequence lifetime of a B2 star - 
the subtype in which the Be star frequency reach its maximum - is 
26 Myr. The lack of Be stars in the youngest clusters indicates that 
a Be star cannot be a young object, but an object close to the end of 
its life in the main sequence. 

However, the above results were drawn from inhomogeneous 
sets of data from the literature. In particular, published ages are of- 
ten affected by important uncertainties. The accurate and homoge- 
neous dating of young open clusters with different Be star content 
is hence a key issue to check the evolutionary hypothesis for the Be 
phenomenon. 

The obtained age and associated error safely exclude an age 
lower than 10 Myr, and places NGC 663 in the age interval (13- 
25 Myr) where the maximum Be star frequency occurs. A similar 
result was obtained in CF02 for the clusters h and x Persei, which 
also have a high Be star content. On the other hand, in a recent pa- 
per. fFuetal]j2003h failed to find any Be star in the younger cluster 
Trumpler 24. The results of this work, together with the two re- 
ferred to above, are pieces of evidence supporting the view that the 
Be star phenomenon is an evolutionary effect appearing at the sec- 
ond half of the main sequence lifetime of a B star. 



6 CONCLUSIONS 

We have presented CCD uvbyP photometry for stars in the central 
area of NGC 663. We have obtained the cluster astrophysical pa- 
rameters from the analysis of the B type star range. The reddening 
is highly variable, with values ranging from E(b-y) = 0.639±0.032 
in the central part to E(b-y) = 0.555 ± 0.038 in the south-east. The 
distance modulus is found to be 11.6±0. 1 mag. (2.1 Kpc), and the 
age log f = 7.4 ± 0. 1 years (25+J Myr). 

A detailed analysis of the Be star population of NGC 663, 
the distribution of Be stars along the main sequence and their fre- 
quency as a function of the spectral subtype will be presented in a 
forthcoming paper. 



ACKNOWLEDGMENTS 

We are grateful to the Observatorio Astronomico Nacional for 
the allocation of observing time in the 1.5m. telescope, and for 
support during observations. This research has made use of the 
WEBDA database, developed and maintained by J.C. Mermilliod, 
the SIMBAD database, operated at CDS, Strasbourg, France, and 
the NASA Astrophysics Data System Abstract Service. This work 
has been partially supported by the Plan Nacional de Investigacion 
Cientifica, Desarrollo e Innovacion Tecnoldgica del Ministerio de 
Ciencia y Tecnologia and FEDER, through contract AYA2000- 
1581-C02-01, and by Generalitat Valenciana AVCyT through grant 
GRUPOS03/170. 



REFERENCES 

Capilla G., Fabregat J., 2002, A&A 394, 479 (CF02) 
Capilla G., Fabregat J., Baines D., 2000, in Smith M.A., Henrichs 
H.F., Fabregat J., eds, Proc. lAU Coll. 175, The Be phenomenon 
in early-type stars, Astron. Soc. Pac. Conference Series 214, p. 
63 

Crawford D.L., 1978, AJ 83, 48 

Fabregat J., Torrejon J.M., 2000, A&A 357, 451 

Fabregat J., Torrejon J.M., Reig P., Bernabeu G., Busquets J., 

Marco A., Reglero V., 1996, A&AS 1 19, 271 
Figueras J., Torra J., Jordi C, 1991, A&AS 87, 319 
Fu J.N., Sterken C, Duerbeck H.W., Mennickent R.E., 2003, 

A&A 412, 97 

Hoag A.A., Johnson H.L., Iriarte B., Mitchell R.I., Hallam K.L., 
Sharpless S., 1961, Publ. USNO, Second Ser., 17, 343 

H0g E., Fabricius C, Makarov V.V., et al. 2000, A&A 355, L27 

McCuskey S.W., Houk N., 1964, AJ 69, 412 

Mermilliod J.C, 2000, in Rebolo R., Zapatero-Osorio R.M., eds. 
Very Low-Mass Stars and Brown Dwarfs, Cambridge Univ. 
Press, Cambridge, p. 3 

Phelps R.L., Janes K.A., 1994, ApJS 90, 31 

Perry C.L., Olsen E.H., Crawford D.L., 1987, PASP 99, 1184 

Pietrzyhski G., 1997, AcA 47, 21 1 

Pigulski A., Kopacki G., Kolaczkowski Z., 2001a, A&A 376, 144 
Pigulski A., Kopacki G., Kolaczkowski Z., 2001b, AcA 51, 159 
Sanduleak N., 1979, AJ 84, 1319 
Sanduleak N., 1990, AJ 100, 1239 

Schaller G., Schaerer D., Meynet G., Maeder A., 1992, A&A 96, 
269 

Stromgrem B., 1966, ARA&A 4, 433 

Tapia M., Costero R., Echevarn'a J., Roth M., 1991, MNRAS 253, 
649 

Torrejon J.M., 1997, PhD Thesis, Univ. Valencia 
van den Bergh S., de Roux J., 1978, AJ 83, 1075 
Wallace PT, 1998, Starlink User Note 5.17, Rutherford Appleton 
Laboratory 

Wallenquist A., 1929, Medd. Astr. Obs. Uppsala n. 42 



© 2004 RAS, MNRAS 000,[T]-?? 



